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Abstract 



We discuss thermodynamic properties of dark energy using the formalism of field 
theory at finite temperature. In particular, we apply our formalism to a purely kinetic 
Q" 1 type of k-essence. We show quite generally that the entropy associated with dark 

5^ ' energy is always equal or greater than zero. Hence, contrary to often stated claims, a 

violation of the null energy condition (phantom dark energy) does not necessarily yield 
a negative entropy. In addition, we find that the thermal fluctuations of a k-essence 
field may be represented by a free boson gas with an effective number of degrees of 
■ freedom equal to c~ 

' 

1 Introduction 

In a number of recent papers (HEJEHHEJEIITIIBIEI UHl El D2i EE], thermal properties of 
dark energy have been discussed based on the assumption that the dark energy substance 
is a thermalized ensemble at a certain temperature with an associated thermodynamical 
entropy. It is usually assumed that this temperature is an intrinsic property of dark energy 
(DE) rather than the temperature of the heat bath fixed by surrounding matter. A popular 
model of DE is the so-called k-essence [H] which was originally introduced as a model for 
inflation [15j. The purpose of this paper is to analyze thermal properties of a grand canonical 
system described by a purely kinetic k-essence at nonzero temperature. Our aim is to extend 
and analyze classical solutions to field equations at nonzero temperature. 

Dark energy is usually assumed to be barotropic, i.e., described by an equation of state 
(EOS) in the form p = p(p). Equivalently, one defines a field theory Lagrangian, in which DE 
is described in terms of a classical self-interacting field coupled to gravity (for a recent review, 
see [16]). Then, the EOS may be deduced from the energy- momentum tensor obtained from 
the variational principle. In order to explain an accelerated expansion the DE fluid must 
violate the strong energy condition (13 [18], which requires 3p+p > together with p+p > 0. 
The so-called phantom models of DE have p + p < 0, thus violating even the null energy 
condition. 
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Clearly, from the EOS alone it is not possible to uniquely determine the thermodynamic 
properties of a system. One simple example is the EOS p = p/3 which may describe a 
massless boson gas at T ^ (hence S 7^ 0) but also a massless degenerate Fermi gas at 
T = (hence 5 = 0). A similar situation occurs for any barotropic EOS. 

A consistent grand canonical description of DE involves the thermodynamic equations 
with two variables: the temperature T and chemical potential p.. The chemical potential is 
associated with a conserved particle number N related to the shift symmetry 9 — > 9 + const. 
As we will demonstrate, the resulting thermodynamic equations do not require negative 
entropy even in the phantom regime when the null energy condition is violated. We show 
that if there exist a nontrivial, stable configuration which we call condensate characterized 
by the pressure p c & and the density p c &, then the pd + p c d term in the expression for entropy 
is precisely canceled out by the particle-number term. The only nonvanishing contribution 
to the entropy is due to thermal fluctuations that yield a thermal ensemble similar to a boson 
gas at nonzero temperature. 

We organize the paper as follows: In section [2] we recapitulate the basic hydrodynamics 
of purely kinetic k-essence. In section [3] we introduce the chemical potential corresponding 
to the conserved charge which is related to the shift symmetry. Basic thermodynamics is 
discussed in section HI The grand canonical and canonical partition functions are derived in 
section [5] where we also include a brief comment concerning a general k-essence. Discussion 
and conclusions are given in section El 



2 Purely kinetic k-essence 

Consider the action 

S = J d A x ^/— ~g 
where R is the curvature scalar and 



C = m'W(X) ; X = <r«y>,„ (2) 

is the Lagrangian for the scalar field 9 of dimension m" 1 . The dimensionless function W 
depends only on the dimensionless quantity X. Such theories have been exploited as models 
for inflation and dark matter /energy, for example a purely kinetic k-essence [TH [15J [191 ED] 
or ghost condensate [2T1 [221 [23] . The simplest nontrivial example of purely kinetic k-essence 
is a ghost condensate Lagrangian [21] 

£ gh = m 4 (l-X) 2 , (3) 

which has been also studied in the context of dark matter/energy unification [20]. Another 
example originates from the Dirac-Born-Infeld description of a <i-brane in string theory with 
the scalar Born-Infeld type Lagrangian 



£ DB i = -m*VT^X. (4) 

This Lagrangian was derived from the Nambu-Goto action for a <i-brane moving in the d+2- 
dimensional bulk [2~4"1 [25] (for a simple derivation, see also [26]). It may be easily seen that 
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(jl]) yields the EOS p oc — p _1 of the Chaplygin gas, an exotic fluid which has been suggested 
as a model for unification of dark energy and dark matter J27J 123 [29] . The generalization 
to p oc — p~ a (0 < a < 1), was suggested [30] and shown to derive from the Lagrangian 

£ gen = -m 4 (l-X^)/ 2 «) Q/(1+Q) (5) 

which represents yet another example of purely kinetic k-essence. Subsequently, the term 
"quartessence" was invented [31] to denote unified dark matter /energy models. 
From (jSJ) the equation of motion for 9 follows 

{Cxg^e,^ = o, (6) 

where Cx denotes the partial derivative with respect to X. Equation (jSJ) implies the existence 
of a conserved current 

f = -C x g^6, v . (7) 
m 

related to the invariance under the constant shift 9 — > 9 + const of the scalar field 9. The shift 
symmetry reflects a correspondence between a purely kinetic k-essence and a U(l) symmetric 
complex scalar field theory in the so-called Thomas- Fermi approximation J28J [32] . It may 
be shown [33] that the Lagrangian (j2j) is equivalent to 

C = rig^®*^^ -m 4 U(ri\$\ 2 /m 2 ), (8) 

if the amplitude of the complex scalar field $ = |$|e im£ ' varies sufficiently slowly. Here rj = 1 
for a canonical scalar field and rj = — 1 for a phantom. The potential U(Y) is related to 
W(X) in (j2J) by a Legendre transformation 



with X and F satisfying 



+ C/(F) = XY, (9) 

Then the current (jTj) corresponds to the Klein-Gordon current 

jUe^^-i^), (12) 

and hence, the conserved quantity N corresponds to the usual U(l) charge of the complex 
scalar field. 

Assuming X > the field 9 may be regarded as a velocity potential for the fluid 4- velocity 



u" = g^9, u /VX (13) 

satisfying the normalization condition = 1. As a consequence, the energy- momentum 
tensor 

= ILxBJr ~ Cg»v (14) 
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derived from the Lagrangian £ in takes the perfect fluid form, 



T/mu = (p + p)u^u v - P9va> , 
with the parametric equation of state 

p = £, (16) 
p = 2X£ x -£, (17) 

and the speed of sound 

4 = tJL - (18) 

Px 

As before, the subscript X denotes the partial derivative with respect to X. A perfect fluid 
description applies only for X > 0. Furthermore, equations (TLB"]) and ffTTj) imply that the 
domains where £x > correspond to a canonical scalar field (p + p > 0) and those where 
£x < to a phantom. In particular, if in the neighborhood of X = 0, £ ~ rjX, the kinetic 
term is canonical for rj = 1 and is of phantom type for rj = — 1. The field that behaves as 
phantom near X = is also known under the name ghost. 

The conserved particle number associated with the current (J7J) is 

N = yVtffi„ = J nu^, (19) 

where the integration goes over an arbitrary spacelike hyper-surface E that contains the 
"particles" . Using the definition (|T3l) for the velocity, we obtain the particle-number density 

n = —VX£ X ■ (20) 
m 

A similar expression, with the right-hand side differing only in a dimensionfull constant 
factor, has been derived also in |12j . 



3 Effective Lagrangian 

In this section we construct the effective Lagrangian as a function of the chemical potential 
fi associated with the conserved particle- number f[T9"j) . In the Hamiltonian formulation [51] 
we choose the hyper surface E at constant time so that the total number of particles (TT91) 
becomes a volume integral 

2 r 1 r 8C 

N — — I £ x g 0u e,„dv = - / ^-dV. (21) 
Next, we define the grand canonical partition function 

Z = Tr e-^-^ = [ [dn] [ [d9] exp f dr [ dV ( m^- - H(tt, 9 A + ^n) , (22) 

J ./periodic JO J \ OT 171 J 

where r denotes Euclidean time. The Lorentzian and Euclidean times are related by 

t = -ir, (23) 
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as usual. In equation fl22|) we employed fl2Tl) to express the particle number N in terms of 
the conjugate momentum field 

' = Bii- (24) 
The Hamiltonian density H, is defined as usual by the Legendre transformation 

W(7r J i )=7r0, o -£(0,o,0,i), (25) 

with 

= (26) 

A formal functional integration of (|22|) over 7r yields the partition function expressed in terms 
of the effective Euclidean Lagrangian 

Z = f [d9] exp - f dr \ dVC E (9, y) . (27) 

J periodic «/ J 

Since the exact functional integration over tt is possible only if 7i is at most a quadratic 
function of tt, we apply the saddle point approximation. In this approximation, the path 
integral is just the integrand evaluated for the field tt which solves the saddle point condition 

i»f* + ii = 0. (28) 
or air m 



Using ff26l) we find 

and with (I25p we obtain 



«, = *^ + ^, (29) 



- C E = m— - 7r9, + + C(6 fl , 6,) = C{i— + ^). (30) 
or m or m 

Hence, the effective Euclidean Lagrangian is obtained from the Lagrangian (T5]) by replacing 
the derivatives of the field 9 by 

e, v ^e, v + ^s° v . (31) 

m 

Note the difference and similarity with the usual Euclidean field theory prescription [35] 

for a canonical complex scalar field $ where the + or — sign is taken when the derivative 
acts on $* or <3>, respectively. 
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4 K-essence thermodynamics 



In a grand canonical ensemble the thermodynamic quantities such as pressure p, energy 
density p and entropy density s are functions of the temperature T and chemical potentials 
Pi associated with conserved particle numbers iVj. For simplicity, we consider only the case 
of a single conserved particle number N and its associated chemical potential p. 
We start from the standard thermodynamical equation 

sT = p + p — pn. (33) 

Clearly, if p = the positivity of entropy requires p + p > 0. Hence, ignoring p one could 
conclude that a phantom field must necessarily yield a fluid with negative entropy. However, 
this conclusion is incorrect since generally p ^ and the entropy density given by (|33|) need 
not be negative. In fact, as we will shortly demonstrate, given arbitrary temperature T it is 
always possible to find a range of p such that s > 0. 

The entropy and particle-number densities may be expressed as partial derivatives of p 



dp 
df 



dp 
dp 



(34) 



Using this and (1331) we also find 



dp dp 

V + P-T- + p-. (35) 



We now apply these general thermodynamic considerations to purely kinetic k-essence. 
Equations (TTB]) and (1171) define a parametric EOS which is essentially barotropic. In other 
words, by eliminating parametric dependence, the EOS may be put in the form p = p(p). Ob- 
viously, with a barotropic EOS alone one can uniquely determine neither T nor p. However, 
the k-essence relation 

p + p = 2XC X} (36) 

which follows from ([TBI) and (fTTj) . may be used to reduce the arbitrariness in functional 
dependence on T and p. From (l3l)|) combined with (13"6T) it follows that the variable X as a 
function of T and p satisfies a partial differential equation 

The most general solution to this equation is a homogeneous function of 2nd degree which 
may be written as 

X = ^f(T/p). (38) 

Here / is an arbitrary positive dimensionless function of x = T / p. However, the consistency 
with (1201) places further restrictions on /. 

The entropy and particle-number densities may be calculated from p using (|T6|) and 
With help of (EHD we find 

8 = ±f'£x; n=(^f-— 2 f')£ x . (39) 
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Combining the second equation with ([20]) we obtain a simple differential equation for / 

xf'-2f + 2^/f = 0, (40) 

with the solution 

f = (Cx + l) 2 , (41) 

where C is a constant. Now we require S = at T = 0, which implies /'(0) = and hence 
C — 0. This in turn implies /(x) = 1, 

*=4> ( 42 ) 

and s = 0. Hence, we come to the conclusion that the thermodynamic quantities, such 
as pressure and density, derived from the classical kinetic k-essence field theory are not 
temperature dependent and can only depend on the chemical potential fi. Besides, we find 
that the corresponding entropy is zero. This result is not unexpected since it is well known 
that classical solutions to a scalar field theory (e.g., nontopological solitons) correspond to 
saddle point solutions of the Euclidean path integral at nonzero /i and zero temperature 
[36], 137]. The existence of stable, nontopological solitons (Q-balls) that have a nonzero value 
of the conserved charge was proven [38J for a complex canonical scalar field theory of the type 
dHJ) for a class of "acceptable" potentials U. Hence, based on the Thomas-Fermi equivalence 
mentioned in section [2J we infer the existence of similar stable classical configurations in a 
corresponding class of purely kinetic k-essence. Otherwise, if such configurations do not exist, 
then the only stable configuration at zero temperature is trivial, i.e., with p = p = n = 0, in 
which case the entropy S is also zero. 

However, if DE interacts with thermalized particles from the surroundings, its thermal 
fluctuations about the classical solutions would be in equilibrium with the heat bath at 
nonzero temperature. In this case we will have S > as usual. In the next section we derive 
the corresponding thermal contribution to the partition function. 



5 Thermal fluctuations 

In this section we derive the grand canonical thermodynamic potential and canonical free 
energy for a self-gravitating kinetic k-essence fluid contained in a sphere of large radius in 
equilibrium at nonzero temperature T = 1/(3. Although the equilibrium assumption implies 
static metric, our analysis may be also applied to an expanding cosmology in which case 
the thermodynamic parameters, such as temperature and chemical potential, are functions 
of the scale parameter a. In this case the background metric g M „ generated by the mass 
distribution is assumed to be a slowly varying function of time on the inverse temperature 
scale. More precisely, we assume 

9 ~W«T- (43) 

We neglect the influence of matter and radiation assuming that their interaction with DE is 
small and serves only to provide a heat bath at the temperature T. 
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5.1 Grand canonical ensemble 



In a grand canonical ensemble we introduce the chemical potential fi associated with the 
conserved particle number N as in section [3j The partition function is given by 

Z = Tr e'^-^ = J \dg] [d6]e- s ^ s \ (44) 

with the Euclidean actions S g and Sk for the gravitational and the k-essence fields, respec- 
tively. The gravitational part may be put in the form (39J, HQ] 



S z = "77- / d^x^R - [ d 3 xVh(K - K Q ), (45) 

107T JY 07T JdY 

where h is the determinant of the induced metric on the boundary, and K — Kq is the 
difference in the trace of the second fundamental form of the boundary dY in the metric 
Qe and the flat metric. The boundary is a timelike tube which is periodically identified in 
the imaginary time direction with period j3. Thus, the functional integration assumes the 
periodicity in imaginary time and the asymptotic flatness of the metric fields. The k-essence 
Euclidean action is 

S k = J d^x^fg^C-E = — J d A Xy/^g C, (46) 

with the Euclidean Lagrangian Ce given by ( |30l) . The path integral is taken over asymptot- 
ically vanishing fields which are periodic in imaginary time r with period (3. The dominant 
contribution to the path integral comes from metric and k-essence fields, which are near the 
classical fields. The metric is assumed to be static spherically symmetric, and asymp- 
totically flat, and is positive definite with the Euclidean signature due to the substitution 

(EHD. 

In order to extract the classical contribution, we decompose 9 as 

d(x) = Q{x)+m-\(x), (47) 

where is a solution to the classical equation of motion which we will call condensate and 
(p describes quantum and thermal fluctuations around 0. The factor m~ 2 is introduced for 
convenience so that the field if has the canonical dimension one. The action Sk splits in two 
parts = S c i[Q] +Sth[ ( f] and if we neglect quantum fluctuations of the metric, the partition 
function factorizes as 

Z = Z g Z c \Z th , (48) 

where Z g = e~ Ss and Z c \ = e~ Scl represent the saddle-point gravitational and classical 
contributions, respectively. The last factor is the contribution due to thermal fluctuations 

Z th = J[d<p)e- S * M . (49) 

First, we calculate the classical contribution. The classical part of the action is given by 

S c i = - f d 4 x^g~£{X). (50) 

Here and below 

X = <T0„0,,, (51) 
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where 

„ .dQ u <96 a 
or m at m 

*<-£■ « 53 > 

By making use of the field equation ([6]) in the comoving frame, i.e., in the reference frame 
where the 4-velocity components are 



it follows 

9 = + i? = -iwr + i? , (55) 

where u and # are constants. Furthermore, using ( 15T1) with ( 1521 and comparing with ( H2l 
we conclude that u = 0. Thus, the background field 6 is constant. However, since 6 is 
given by (|52|) . the quantity X may still be a function of x and a slowly varying function of 
time. Equations (|5T|) - (!52|) give 



,2 r.2 



X = g 00 ^ = ^ , (56) 

where /2 = /i/y^oo is the local chemical potential. This equation is consistent with equation 
( H2"j) obtained from purely thermodynamic considerations. 

In a spherically symmetric static geometry, regular solutions to equation ([6]) coupled 
with Einstein field equations describe dark energy stars [HJ H2] at zero temperature. A 
particular example of Born-Infeld type k-essence stars have recently been studied in jl3]. In 
the context of cosmology, the background is spatially homogeneous isotropic configuration 
with X depending on time through cosmological scale dependence. 

According to our assumption ( fl3l the variation of X is small on the time scale comparable 
with the inverse temperature, i.e., dX/dt -C Xj ' (3. Then, the condensate contribution to the 
partition function may be written as 

lnZd = (3 J^d 3 x^g~ £(X) (57) 

where £ is a spacelike hyper-surface that contains the condensate. Using this equation we 
find the net number of particles in the condensate 

1 dlnZ cl f 3 2/7 

N = -— = = / a x^/g^ -^C x (58) 



j3 d/j, Jt, m 2 
where g^ = det(— gry), i,j = 1,2,3. Alternatively, the covariant definition ([19]) yields 

N= f d 3 x^/g^n, (59) 



where n is the particle-number density in the condensate. Therefore, we identify the particle- 
number density due to the condensate as 

n = 2-^£ x (60) 
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which coincides with f )20|) as it should. 

Next, we calculate the thermal contribution to Z starting from equation (j4"9j) where the 
action Sth[</?] is derived from C(X) by keeping only the quadratic term in the expansion of 
C in powers of <p. We find 



Sth = -m" 4 J fay/^rd&dvtp , (61) 

where 



r = £x<T + 2£ xx g» a g^e, a e, (3 . (62) 

The symbol <9 M in equation f[6"Tj) denotes the partial derivative with respect to a Lorentzian 
coordinate x^. In particular, <9 = d/dt = id/ dr. The /i dependence is absorbed in f^ u 
through the prescription (|52|) . The action describes a massless scalar propagating in an 
effective (or emergent) acoustic geometry [HI H5j HHJ H7J provided the field equation for ip 

= (63) 

is hyperbolic, i.e., provided the effective metric tensor f^ u has the Lorentzian signature. This 
holds if and only if 

detr = C x g- 1 c- 2 (64) 

where 

_ 2 C x + 2£xx 

c s = r ^ 65 ) 

l~,x 

is the inverse speed of sound squared. It may be easily verified that the same equation is 
obtained from the standard hydrodynamic definition 

2 d P 

c. - 



(66) 

s/n 



dp 

Since g < 0, is Lorentzian if and only if c 2 s > 0. Hence, the hyperbolicity condition ([6] 
is equivalent to the requirement of hydrodynamic stability. 

In addition, it seems physically reasonable to require c 2 s < 1 in order to avoid possible 
problems with causality. However, it has been shown [49l [50] that if k-essence is to solve 
the coincidence problem there must be an epoch when perturbations in the k-essence field 
propagate faster than light. In cosmology and astrophysics, it is usually assumed on the basis 
of causality that the speed of sound cannot exceed the speed of light [171 EH EE] . Hence, it 
has been argued [19] that k-essence models which solve the coincidence problem are ruled 
out as realistic physical candidates for dark energy. In contrast to this, it has been argued 
[UJ EH E2] that superluminal sound speed propagation in generic k-essence models does 
not necessarily lead to causality violation and hence, in spite of the presence of superluminal 
signals on nontrivial backgrounds, the k-essence theories are not less legitimate than General 
Relativity. 

As we are not presently interested in discussing or solving the coincidence problem we 
stick to < c 2 s < 1. This condition is fulfilled if both C x and C xx are simultaneously 
either positive or negative. In the latter case from (1361) we have p + p < and hence, 
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our thermodynamic analysis is not restricted only to models that satisfy the null energy 
condition. 

Assuming < <? s < 1, equation ( joTT) may be put in the form 

5th = ~\l <*W=G G^d^ip . (67) 
2 Jy 



where 

c Tn 1 

= 2 C -^—[g^ - {I - -)uV] . (68) 
Cx d 



The matrix G^ u is the inverse of the acoustic metric tensor defined as 



2 c s m 

where 

C- 

u„ = - 

'X 



1 c 

G »v = 7;—^[9»v-(l-c 2 s)u ll u v \, (69) 



= % ; = (70) 



is the velocity of the fluid. 

The determinant G is given by 

16£ 4 

G = detG, u = -^g. (71) 

Since the metric is static by assumption , i.e., g^ u is independent of t and goi = 0, the same 
is true for the acoustic metric G,„, and the determinant factorizes as G = —GnnGt^ where 



Cr(3) = — det Gij] i,j = 1, 2, 3. By making use of the substitution r = ryGoo we obtain 

d 3 x^G {3) / df£(f,x), (72) 



where 

1 



C(f,x) = -\{d fV f - G ij diipdjip] . (73) 



Here we have introduced the effective inverse temperature 



P =y /cW = J^P, (74) 

where the quantity /3 = y/goo/3 is the usual local inverse temperature. Equation ( 17"4|) is the 
well-known Tolman condition for thermal equilibrium in curved space [531118]. The parameter 
(3 retains its usual interpretation as the asymptotic value of the inverse temperature [51] . 

Using the approach developed in [37] the path integral in (149]) may be easily calculated. 
We obtain the thermal part of the partition function in the form 



lnZ th = - jj^x^G^j A^ln(l-e^). (75) 
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This expression may be regarded as a proper volume integral 



In Z th = - J^d 3 x y/g^j — In z (76) 



of the local partition function 
\n.z = -V 



\ 9(3) J (2tt)3 v ' 360 \ (3) 



G (3) 3-3 



p, (77) 



from which the pressure, energy density, and entropy density may be derived in the usual 
way: 

Pth = In z = — c7 3 p-\ (78) 
y VP 360 s P ' V ; 

Pth = -I-lln2 = — c7 3 ^- 4 , (79) 
m Vdp 120 s P ' V ; 

^ = /5(pth + Pth) = ^c; 3 ^- 3 . (80) 

These expressions are the well-known thermodynamic equations that represent a massless 
relativistic Bose gas in curved space. The factor c~ 3 is the effective number of degrees of 
freedom which depends on the chemical potential for a particular k-essence model. This 
shows that a sizable fraction of radiation today can be attributed to thermal fluctuations of 
DE. For example, for the Born-Infeld type k-essence (j3j) we find 

cf = (l-S/mr 3/2 . (81) 

If we use this model for a dark matter/energy unification, the fit with homogeneous cosmol- 
ogy today would yield the speed of sound squared of the order c 2 s = Qa with the effective 
number of degrees of freedom f2 -3 / 2 w 1.6, comparable with 2 for photons. 

The gravitational part of the partition function may be calculated from (1451) with help 
of Einstein field equations. Using the result of Gibbons and Hawking for the surface term 
[391. we obtain 



\nZ g = -[3M+ / rf 4 Xv^T °, (82) 

where M is the total mass and is the energy-momentum tensor of k-essence field averaged 
with respect to the partition function fj48|) . The averaged energy-momentum tensor may be 
split up into two parts: 

Tfju, = T c d [if 4" 7th \ivi (83) 

where the first term on the right-hand side is the classical part which comes from the con- 
densate and the second term represents the thermal fluctuations. 

It may be shown that both terms are of the form (fl5"|) which characterizes a perfect fluid. 
Hence, we have 

Tq = ped + Pth ; T- = p cd + pth , (84) 

where the thermal pressure and energy density are given by (ITS"]) and f!79|) . respectively. The 
condensate pressure p c d and energy density p C( j are given by fTTBT) and ffTT|) . respectively, in 
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which the quantity X is given by fl56|) . Putting the condensate floTj) . the thermal (1761) . and the 
gravitational fl82l) contributions together, we find the total grand canonical thermodynamic 
potential as 



Q(J3, At) = -- InZ = M - f d 3 



%V^g (pcd + Pcd + pth + Pth) • 



(85) 



This is just the standard form of the thermodynamic potential 



n = E-T3-fiN, 



(86) 



in which we identify the energy E with the total mass M, the TS term with the thermal 
contribution 



and the /j,N term with the contribution of the condensate 




Note that the entropy S is strictly positive because the right-hand side of ( 1H7|) is positive 
and the quantity T is a positive temperature of the heat bath. 

5.2 Canonical ensemble 

Now consider a self-gravitating k-essence fluid with N particles contained in a two-dimen- 
sional sphere of large radius in equilibrium at nonzero temperature T = 1/(3. In a canonical 
ensemble, instead of the chemical potential \i we fix the particle number. Thus, a canonical 
ensemble is subject to the constraint ( |59l) with iV fixed. 

The free energy of a canonical ensemble may be derived from the grand canonical partition 
function with the help of the Legendre transform 



The quantity jj, in this expression is an implicit function of N and T, such that for given iV 
and T the constraint (J5l?j) is satisfied [55J. From flHBI) and (J55|) with fl6"Ul) it follows that 



The second term on the right-hand side is related to the entropy through ( !80l) and the free 
energy may be expressed in the familiar form 




(87) 



F((3,N)=Q((3, f i)+ f iN. 



(89) 




(90) 



F = M -TS, 



(91) 



where the total entropy S is defined as a proper volume integral 




(92) 



of the entropy density s given by (J8*0|) . 
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5.3 More general k-essence 

The above considerations may be similarly applied to a general class of k-essence models 
described by 

c(e,x) , 



S = Id x \J—g 



(93) 



16nG 

with the most general Lagrangian which, in addition to X, depends explicitly on the scalar 
field 9. The hydrodynamic quantities u^, p, p, T^ u , and c s associated with ([93"]) are, as before, 
defined by (I13p - (|18p . However, unlike in purely kinetic k-essence, the equation of motion for 
9 

dC 

QCxgrOr)* ~ W = (94) 

is no longer of the current conservation form and the analysis of sections [3] and H] does not 
apply Here, for example, we cannot establish a correspondence with the canonical complex 
scalar field theory. As in the case of a real scalar field, because of the absence of a conserved 
Noether current there exist no stable nontopological solutions [56j although time dependent 
solutions similar to oscillatons [57J [58] , or static unstable configurations similar to unstable 
scalar solitons [59], are not excluded. The only stable solutions at zero temperature are 
trivial, i.e., those with p = p = so that the condensate contribution to the partition 
function is absent. 

As there is no conserved particle number the canonical and grand canonical ensembles 
coincide, i.e., F(/3) = fi(/3), and the free energy F is given by (1901) . The thermal contribution 
is of the form (j4"9l where the action <Sth[<^] is obtained by expanding C(X, 9) in powers of cp 
with 9 = + m~ 2 ip, as before. Here, describes a configuration for which p = p = 0. Using 
the procedure described in section 15.11 we find the expression similar to ([6j 



Sth = ~\ Iy d4x ^ {GTd^ip - m 2 ^ 2 ) , (95) 

where 

m 2 cS = ^| (2X£ xee - C @e + 2^6^ (96) 

is the effective mass jU], X = g^Q^Q^, and the matrices f^ v and G^ u are defined as in 
f [62|) and (168]) . respectively. The subscripts X and in (196]) denote the partial derivatives 
with respect to X and 0, respectively. The quantity c s is the sound speed defined as in ( 165]) . 
The action ( 195]) describes a massive scalar propagating in an effective acoustic geometry 
provided the condition (1641) is met. Again, using the approach developed in [37] we find the 
expression for the thermal partition function in the form of the proper volume integral (176]) 
over the local partition function 



ln, = -y ^/^ln(l-e^) (97) 
\ 9(3) J (2tt) 3 

where E = yj + q 2 and g^) and are the determinants of the respective spatial metrics 
defined as in section 15.11 From this partition function we find the usual expressions for the 
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pressure, energy density, and entropy density of an ideal gas of massive bosons 

9eS f d 3 q „-pE\ 



/d 3 q 1 
( 27 r)3 1 - e -!>E > (") 

s = PiPth + Pth) (100) 
with the effective number of degrees of freedom g e g — l/c% which depends on the model. 



6 Conclusions 

We have derived a grand canonical and canonical description of k-essence type of DE. The 
thermodynamic equations are generally expressed in terms of two variables: the temperature 
T and chemical potential \i. The chemical potential is associated with a conserved particle 
number N related to the shift symmetry. The derived thermodynamic equations show that 
the entropy does not have to be negative even in the phantom regime, contrary to the claims 
often stated in the recent literature (see, e.g., [8] and references therein) that a violation of 
the null energy condition implies negative entropy. We have demonstrated that the entropy 
is greater or equal to zero and is strictly zero at zero temperature. We have show that if 
there exist a nontrivial, stable configuration which we call condensate characterized by the 
pressure p c d and the density p C( j, then the particle- number term in the expression for entropy 
cancels out the contribution of p c ^ and p cc j. Furthermore, we have shown that the only 
nonvanishing contribution to the entropy is due to thermal fluctuation analogous to those 
of a massless boson field. The thermal ensemble behaves as a free massless gas at nonzero 
temperature with an effective number of degrees of freedom equal to c~ 3 . Similarly, thermal 
fluctuations of a general k-essence field yield an effective free gas of massive bosons. 
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